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Abstract — Bliss schemes of a run length limited (RLL) codec 
in combination with an LDPC codec, generate LDPC parity 
bits over a systematic sequence of RLL channel bits that are 
inherently redundant as they satisfy e.g. a d — 1 minimum run 
length constraint. That is the subsequences consisting of runs 
of length d — 1, viz. . . . 010 . . . and . . . 101 . . ., cannot occur. 
We propose to use this redundancy during LDPC decoding in a 
Bliss scheme by introducing additional d-constraint nodes in the 
factor graph used by the LDPC decoder. The messages sent from 
these new nodes to the variable or codeword bit nodes exert a 
"force" on the resulting soft-bit vector coming out of the LDPC 
decoding that give it a tendency to comply with the d-constraints. 
This way, we can significantly reduce the probability of decoding 



Index Terms — Bliss schemes, RLL codes, LDPC codes, factor 
graph, modified concatenation. 



I. Introduction 

BLISS schemes [1], also called modified concatenation 
schemes by Fan [2], place a modulation encoder and 
decoder as outer parenthesis and a systematic ECC encoder 
and decoder as inner parenthesis around the storage channel. 
This way, the modulation decoder follows the ECC decoder, 
rather than precedes it, thus, avoiding error propagation by the 
modulation decoder at the input of the ECC decoder as in the 
standard concatenation scheme [2]. As shown in Fig. 1, the 
parity bits generated from the sequence of modulation encoded 
channel bits are encoded with a second modulation encoder, 
where both modulation encoders are of the run length limited 
(RLL) type [3] . 

The so called IT '-precoder in Fig. 1 performs an integration 
modulo-2 in order to convert differential bits, also called 
d, k-bits within the context of RLL coding, to so called 
unipolar bits, that are indicative of the type of run on the 
storage medium. The channel SISO detector is shared between 
the systematic and the parity part of the LDPC codeword 
and produces soft-decision values for the unipolar bits. The 
equalizer filter shortens the channel impulse response. The 
SISO RLL decoder operates on these soft-decision unipolar 
bits produced by the channel SISO detectors. Hence, the 
number of states of the RLL SISO decoder is twice the number 
expected if differential bits were used, as the state incorporates 
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an additional polarity bit 1 . 

Bliss schemes have pros and cons in comparison with a 
standard concatenation scheme: 

• Pro: In the Bliss scheme, the bit error rate (BER) at the 
input of the LDPC decoder, at least for the dominating 
systematic part, is the bare BER after bit-detection, which 
has not been multiplied by the error propagation factor 
of the RLL SISO decoder. 

• Con: The error correction capability of the LDPC code is 
characterized by its code rate Rldpc ( ar >d secondarily by 
its block length N). The RLL encoding of the systematic 
part increases the number of bits into the ECC parity 
generator. For instance, for the same LDPC code rate 
Rldpc with a rate 2/3 RLL code, as applies for the 
d = 1 case, the parity generator produces 1.5 times more 
parity bits, due to the RLL encoding of the systematic 
bit sequence. This additional amount of parity lowers the 
net rate of the Bliss scheme. 

In order to mitigate the rate loss of Bliss schemes due to the 
parity generation over redundant data as mentioned above, Im- 
mink [7] introduced a compression codec in Bliss' scheme. For 
storage channels with jitter noise on the transitions between 
run lengths, Zhang et al. [8] reduced the error propagation 
of the compression decoder for a d = 5 RLL code by the 
use of Gray labeling. In this contribution, we take a different 
approach. We accept the redundancy of the systematic channel 
bit sequence, and modify the LDPC decoder to exploit this 
redundancy during its decoding iterations, in order to improve 
its error correction performance. 

Our interest is in the use of the same RLL codes by both 
modulation encoders and low density parity check (LDPC) 
codes [4] as ECC codes with a high code rate (e.g. Rldpc = 
0.90 — 0.95). Especially, in near-field optical storage as en- 
visaged for the fourth generation of optical recording, RLL 
constraints with a minimum run length constraint of d = 1 
are still quite popular as evidenced by the recent RLL code 
design of Coene et al. [5]. In [6], we discussed a stitching 
technique to connect a systematic part of a codeword to its 
adjacent parity part. 



'The use of a unipolar RLL SISO as in Fig. 1 gives a moderate bit error 
performance improvement over a channel SISO that outputs soft differential 
bits and an RLL SISO that with half the number of states. Even better bit error 
performance is possible by combining the channel SISO and the RLL SISO 
in one joint SISO with the product space of the state space of the channel 
SISO and RLL SISO. 
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II. Min-Sum LDPC Decoding Revisited 

The min-sum LDPC decoding algorithm [9] is a simplifi- 
cation of the sum-product algorithm that uses minimum and 
summation operations instead of multiplication and summation 
operations. Both algorithms are a special case of the message 
passing algorithm [10] that passes messages along the edges 
of a so called factor graph. These messages convey soft- 
decision information, e.g. in the form of log-likelihood ratios 
(LLR). The min-sum algorithm always uses a log-likelihood 
representation of the messages. For a binary random variable 
with probability p (pi) of taking on the value (1), the log- 
likelihood equals log(p /pi)- With a properly chosen scaling 
factor after the minimum operation, the performance loss of 
the min-sum algorithm w.r.t. the sum-product algorithm is 
minimal [11] [12]. The choice of the codeword bit node or 
"variable node" degrees in the factor graph [10] of the LDPC 
code seems to have only a minor influence for high LDPC code 
rate. 

Let N be the LDPC 2 codeword length and M be the number 
of parity check equations. The parity check matrix H consists 
of M rows and N columns with elements from the binary 
Galois field GF{2). For a parity check equation with index 
to, < to < M, define the set M(m) of codeword bit 
positions that it checks, i.e. 

M{m) = {n | n = 0, 1, . . . , N - 1, H mn ± 0). 

Similarly, for a codeword bit position n, < n < N, define 
the set M(n) of indices of parity check equations that check 
the given bit position n, i.e. 

M(n) = {to I to = 0, 1, . . . , M - 1, H mn ^ 0). 

The factor graph associated with the parity check matrix H 
has as set of vertices V, the union of the set of N bit nodes 
V var and a set of M parity check nodes V c heck- The set of 
edges E consists of all edges (to, n) for which H mn ^ 0. For 
the sake of ease of discussion, from hereon, we assume that 
the LDPC code is regular, which means that all sets M{m) 
have the same size K, and that all sets M{n) have the same 
size J. 

In general, the maximum likelihood or hard-decision es- 
timate of a binary random variable is determined by the 
sign of its log-likelihood value. The absolute value of a log- 
likelihood value is a reliability of the hard-decision estimate. 
A log-likelihood value of zero corresponds to an erasure (no 
information). 

We now discuss the various operators involved in the min- 
sum algorithm. The operator following VAR that is applied 
inside the codeword bit nodes or "variable nodes" combines 
a number of sources of information (log-likelihood messages) 
about the bit associated with a variable node. These sources of 
information are assumed to be statistically independent. This 
assumption of statistical independence^ translates into a sum 

2 Our method may also apply to other, similar codes with a low density 
parity check matrix, such as repeat-accumulate (RA) codes. 

3 Due to the presence of cycles in the factor graph, this is only an 
approximation. 



operator on the input log-likelihoods Aj, i = 0, 1, . . . , J — 1 
of a variable node. 

VAR(A , Ai, . . . , Aj_i) = A + Ai + . . . + Aj_i. 

The operator CHK that is used inside the check nodes approx- 
imates the log-likelihood value of the exclusive-or "ffi" of its 
presumably statistically independent input variables. 

|CHK(A ,A 1 ,...,A K -i)| 
= amiri(|A |, |Ai|, . . . , \A K -i\) (l) 

signBit(CHK(A ,A 1 ,...,A x _ 1 )) 

(2) 

= signBit(An) © signBit(Ai) © ... © signBit(A X -i) 
where < a < 1 is the aforementioned scaling factor and 

( if x > 0, 
signBit(.x) = | 1 . fx<Q 

For a given iteration of the min-sum algorithm, we define 
the following variables. As usual, in the message passing 
algorithm, messages are sent along edges of the factor graph. 

• t n - The decoder input message into variable node n. 

• u mn - The message sent from variable node n to check 
node to. It is obtained as a function VAR of message t n 
and the last received messages v mn of all check nodes 

to', to' 7^ TO, 

= VAR( t n , (vm'n \ m! G M(n)\{m}) ). (3) 

• v mn - The message sent from check node m to variable 
node n. It is obtained as a function CHK of the last 
received messages u mn of all variable nodes n' , n' ^ n, 

Vmn = CHK( (u mn > | n! S Af(m)\{n}) ). (4) 

• w n - The decoder output messages. Unlike the messages 
u mn the decoder output message w n uses all available 
information in a variable node n. It is obtained as the 
function VAR of the message t n and the last received 
messages v mn of all check nodes to, 

w n = VAR( t n , (v m > n I to' G M-(n)) ). (5) 

A classical implementation of the min-sum algorithm stores 
all received messages. During the first half-iteration, all mes- 
sages u mn are sent from all variable nodes to the check 
nodes. During the second half-iteration, all messages v mn are 
sent from all check nodes to the variable nodes. A constant 
number of iterations can be used. The decoder output messages 
need not be evaluated for all iterations, but only for the final 
iteration. 

III. d = 1-CONSTRAINT NODES IN THE FACTOR GRAPH 

Our approach is to modify the LDPC decoder so that it can 
use the knowledge that in the systematic channel sequence the 
subsequences 01 p and 10 p l, p = 1, 2, . . . , d do not occur. 
Here, a d denotes a sequence consisting of d copies of the 
value a, a = 0, 1. This is achieved by addition of so called 
d-constraint nodes in the LDPC code's factor graph. As stated 
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before, our aim is to improve the error correction capability 
of the LDPC decoder. 

In the definition of a new, modified factor graph for the 
LDPC decoder in a Bliss scheme, it is very convenient that 
the d-constraint is a local constraint [10]. That means, that 
the constraint involves only a fixed, small number d + 2 of 
codeword bits (codeword bit nodes). This sparseness of the 
connectivity matrix between the nodes is essential for efficient 
and effective operation of the modified LDPC decoder. This is 
all the more true, the smaller the value of d, d > is. For this 
reason, and for the practical importance of their applications 
[5], from hereon we concentrate on the d = 1-based Bliss 
schemes. 

The new factor graph with the d-constraint nodes is shown 
in Fig. 3. In general, the degree of a d-constraint node equals 
d + 2, as these nodes need to be able to detect the presence 
of the subsequences 01 d and 10 d l of length d + 2. Define 
the following additional variables in the LDPC decoding 
algorithm. 

• a np - The input message to d-constraint node p sent from 
variable node n. 

• b np - The output message of d-constraint node p sent to 
variable node n. 

Here, the index n of the variable node (LDPC codeword 
variable node) is implicitly restricted to the systematic bit 
positions as there is no d-constraint on the parity part of an 
LDPC codeword. 

As a general principle in message passing algorithms, the 
output message along a certain arc from a certain node 
is only allowed to depend on the (most recently) received 
input messages via all other arcs into that node . Hence, 
for d = 1-constraint node with index n, and most recently 
received log-likelihood messages a„_i !rl , a„,„ and a n +i, n , 
we need to specify three output functions COo,COi and CO2 
in order to generate as many log-likelihood output messages 
b n -i, n , b n ,n, b n+liU , such that 

b n -i,n = COo(a„, n , a n +i, n ) 

b n +l,n = C02(an-l,n, 0,-n,n)- 

Observe that, if one of the input log-likelihood values 
a-n-ii a n, »n+i of a d = 1-constraint node is zero, there 
is no indication that the d = 1-constraint is violated, as it is 
not possible to conclude to a violation from the knowledge 
of fewer than d + 2 = 3 hard-decision estimates. Then, the 
output log-likelihood messages are to be zero. The output 
messages do not need to exert a force on the solution of the 
LDPC decoder in that case. On the contrary, if the input log- 
likelihoods have large absolute values, and their hard-decision 
values indicate a violation of the d-constraint, the output 
values, also should have large absolute values. However, if the 

4 Note, that a d = 1 constraint itself does not imply a non-trivial Hamming 
distance D, D > 1. For instance the d = 1 constrained sequences 
00 . . . 01100 ... 00 and 00 . . . OHIO ... 00 evidently have Hamming dis- 
tance D = 1. However, for the kind of intermediate range of target bit 
error rates after LDPC decoding in our storage application 5 a high minimum 
distance of the LDPC code is not required for good or even superior error 
correction performance, anyway. 



hard-decision values indicate compliance with the d-constraint, 
the output log-likelihoods should be zero. Hence, in the spirit 
of the min-sum algorithm we choose to let 

\CO r {x, y)\ € {0, min(|x|, \y\)}, r = 0,l,2. 

Similar to Eq. (2), one can use a scaling factor to post-multiply 
the above minimum values. 

Also the sign of COo is not allowed to depend on a n _^„, 
etc. The signs are chosen such that they enforce the disap- 
pearance of violations of the d = 1 constraint in the decoded 
LDPC codeword, see Tables I, II and III. As stated before, a 
zero entry in these tables applies when there is no violation 
of the d-constraint. 

IV. Min-Sum LDPC Decoding with 

d = 1-CONSTRAINT NODES IN THE FACTOR GRAPH 

With reference to the above definitions of log-likelihood 
messages, one iteration of the extension of the min-sum LDPC 
decoding algorithm with d = 1-constraint nodes in the factor 
graph is given by the following equations. 
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(6) 



V. Simulation Results 

We experimented with two pseudo-random regular LDPC 
codes. The codeword bit nodes or variable bit nodes all had 
degree J — 3. The shorter LDPC code had a code rate of 
Rldpc = 0.906 and a code length of N = 1728. The bit 
error rates after LDPC decoding obtained from simulations 
using this shorter code are shown in Fig. 4 and 5. In Fig. 4, 
the peak signal-to-noise ratio (PSNR) of white additive noise 
at the channel output is varied for a fixed channel bit length 
of 53 nm, scaled to the numerical aperture and wavelength of 
Blu-ray disc [13]. The user capacity of the disc scaled to the 
physics of Blu-ray disc then equals 30.5 GB. The simulations 
used a channel model of the optical storage readout channel 
using the channel modulation transfer function from the Braat- 
Hopkins formalism [14]. The channel bit length is scaled to the 
physical readout parameters of Blu-ray disc. In Fig. 5, channel 
bit length is varied, at a fixed PSNR of 35 dB. The d = 1 RLL 
code of [5] is used, that has code rate Rrll = 2/3. 
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The number of LDPC iterations was set at 16. We use 
the schedule of Yeo et al. [15], wherein a single check 
node or rf-constraint node update is followed by variable 
updates of the connected variable nodes as a kind of "mini 
decoding iterations." We process the check nodes and the d- 
constraint nodes in sequence. An LDPC decoding iteration 
is then complete when all check nodes and all d-constraint 
nodes have undergone an update. This update schedule about 
halves the number of iterations that need to be performed in 
a decoder. Up to half a million LDPC codeword have been 
simulated per point of the graph. 

In all graphs presented in this section, the top curve shows 
the BER after bit-detection. The middle curve shows the BER 
after standard LDPC decoding. The bottom curve shows the 
BER with decoding using the additional RLL-constraint nodes 
in the factor graph. 

The longer LDPC code has a code rate of Rldpc — 0.955 
and a code length of TV = 6912. The BER simulation results 
for this longer code are shown in Fig. 6 and 7. 

For a physical sector size of the storage medium of approx. 
1 Mbit, the shorter LDPC code length allows combination with 
an 10-bit outer RS code with about maximal length (2 10 — 1 
symbols). The longer LDPC code allows combination with an 
8-bit outer RS code of about maximal length (2 8 — 1). 
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VI. Conclusion 

We conclude that the use of a d = 1 -constraint during the 
decoding of an LDPC -based Bliss scheme, brings an advantage 
in PSNR of around 0.25 dB. Using this technique, the bit 
length of a simulated model of an optical storage channel can 
be decreased by ca. 0.6 %. This suggests that the use of d = 1- 
constraint during LDPC decoding can enable an increase in the 
storage density of around 0.6 %. 



Acknowledgment 

These results were obtained during a joint project between 
Philips Research and Sony. The author is indebted to his 
former colleague Haibin Zhang, now with TNO Telecom, for 
the construction of the LDPC codes used in this research. 
Furthermore, the author would like to acknowledge stimulating 
discussions with his former colleagues Stan Baggen, Wim 
Coene and Bin Yin as well as his project members from Sony, 
Seji Kobajashi, Toshi Horigome, Makato Noda and Hiroyuki 
Yamagischi. 

References 

[1] W. G. Bliss, "Circuitry for performing error correction calculations on 

baseband encoded data to eliminate error propagation," IBM Techn. 

Disci. Bui, vol. 23, pp. 4633-4634, 1981. 
[2] J. L. Fan, Constrained Coding and Soft Iterative Decoding, Norwell, 

MA: Kluwer Academic Publishers, 2001, pp. 144-145. 
[3] K. A. S. Immink, Codes for mass data storage systems, 2nd ed., 

Eindhoven, The Netherlands: Shannon Foundation Publishers, 2004, pp. 

51-64. 

[4] R. G. Gallager, Low-Density Parity-Check Codes, Cambridge, MA: MIT 
Press, 1963. 

[5] W. M. J. Coene, A. Hekstra, B. Yin, H. Yamagishi, M. Noda, A. 
Nakaoki, and T. Horigome, "A new d=l, k=10 soft-decodable RLL 
code with r=2 RMTR-constraint and a 2-to-3 PCWA mapping for DC- 
control," in Proc. of SPIE, Vol 6282, Optical Data Storage 2006, R. 
Katayoma, R. Schlesinga, editors. 




Andries P. Hekstra (M'00) was born in Breda, the 
Netherlands in 1961. He received the "Ingenieur" 
degree in Electrical Engineering from Eindhoven 
University, the Netherlands, summa cum laude in 
1985 with a specialization in multi-user informa- 
tion theory. In 1985-86 he was a Young Graduate 
Trainee at the European Space Agency where he 
worked on spread spectrum telemetry systems in 
Darmstadt, Germany. Successively, he was a Ph.D. 
student at the Electrical Engineering Department of 
Cornell University, Ithaca, USA. There, he studied 
abstractions of VLSI packing problems. In 1990 he joined KPN Research 
in Leidschendam, the Netherlands, where he finished his Ph.D. degree in 
1994 at Eindhoven University of Technology with his professor from Cornell 
University as second promoter. From 1995 to 2000 he investigated automatic 
assessment of video and speech quality using models of human perception 
and cognition. During 2001-2006, he worked at Philips Research, Eindhoven, 
mainly on error correction for optical storage systems. Since Sept. 2006, he 
works on applied communication and information theory at large within the 
research department of NXP Semiconductors, Eindhoven. 



5 



LDPC 
Parily 




RLL 
Encoder 




Precoder 


- 1 


Generator 











&71 — 1,77 ^77+1,77 


bn,n — COi(a n -i in) a n +i t n) 


+ + 


+ 


+ 





+ 






System alic Part 



TABLE I 

Sign of the output log-likelihood message &„,„ of the 
d-constraint node with index i as determined by the function 
COi(a„_i,„, a n +i,ti) 



SISO 
Channel 
Detector 



Fig. 1 . Block diagram of a Bliss scheme. 
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TABLE III 



Sign of the output log-likelihood message b n -x,n of the 
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Fig. 2. Factor graph of an LDPC code. Variable nodes are located in 
the middle. At the left are the input nodes, representing the channel output 
messages. Check nodes are situated at the right. 
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Fig. 3. Factor graph of an LDPC code with additional o!-constraint nodes 
for the case of d = 1. 



Fig. 4. Simulated BER curves for an LDPC-based Bliss scheme with the 
J = 3-regular LDPC code is Rldpc = 0.906 and the code length N = 
1728. 



6 




Fig. 5. Similar graph as Fig. 4, where the peak signal-to-noise ratio (PSNR) 
is fixed at 35 dB, and the bit-length of the channel bits is varied. 




Fig. 6. Similar graph as Fig. 4, with a longer LDPC code (length 6912 bits) 
of a higher code rate 0.955. 
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Fig. 7. Similar graph as Fig. 5, also with a longer LDPC code (length 6912 
bits) of a higher code rate 0.955, as in Fig. 6. 



